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Abstract —Recent work has shown that adaptive CSMA al¬ 
gorithms can achieve throughput optimality. However, these 
adaptive CSMA algorithms assume a rather simplistic model 
for the wireless medium. Specifically, the interference is typically 
modelled hy a conflict graph, and the channels are assumed to 
he static. In this work, we propose a distributed and adaptive 
CSMA algorithm under a more realistic signal-to-interference 
ratio (SIR) based Interference model, with time-varying channels. 
We prove that our algorithm is throughput optimal under this 
generalized model. Further, we augment our proposed algorithm 
by using a parallel update technique. Numerical results show that 
our algorithm outperforms the conflict graph based algorithms, 
in terms of supportable throughput and the rate of convergence 
to steady-state. 

L Introduction 

A central problem in wireless networks is the design of 
efficient link scheduling algorithms in the presence of in¬ 
terference. In the design of scheduling algorithms, there are 
three key performance metrics of interest. The hrst among 
them is the achievable throughput region. The throughput 
performance of a scheduling algorithm is characterized by 
the largest set of arrival rates under which the algorithm can 
stabilize the queues in the network. Secondly, the average 
delay incurred by the packets in the queue should be small. 
The third metric of interest is the computational and commu¬ 
nication complexity involved in implementing the algorithm. 
Scheduling algorithms with low computational complexity and 
low communication overheads are preferable. 

A. Related Work 

A large part of the existing literature on scheduling is 
based on the maximum weight scheduling algorithm ID, 
which is known to be throughput optimal under fairly general 
conditions. However, maximum weight scheduling generally 
requires solving an NP-hard problem during each scheduling 
instant, and is difficult to implement in practice. Eurther, it is 
not directly amenable to a distributed implementation. Several 
low complexity alternatives HI have been proposed but they 
achieve only a fraction of the capacity region, and are hence 
not throughput optimal. 

On the other hand, there are simple random access tech¬ 
niques such as Aloha, CSMA (Carrier Sense Multiple Access) 
which can be implemented in a distributed manner. A dis¬ 
tributed algorithm was developed in Q to adaptively choose 
the CSMA parameters so as to achieve throughput optimality. 
Central to this algorithm is the so called Glauber dynamics, 
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which is a Monte Carlo Markov Chain sampling technique lH, 
0. Specihcally, it is a Gibbs sampler 0 based algorithm. 

A main shortcoming of the existing papers on adaptive 
CSMA is that the results are derived based on rather simplistic 
models for the wireless channels and the interference. Typical 
modelling assumptions used include: 

a. Conflict graph interference model: The interference is 
modelled by a conflict graph or protocol model 0, where 
the transmissions from two links fail, if the links share an 
edge in the conflict graph. In reality however, the success 
or failure of a link depends on the aggregate interference 
from all the active links in the interference range. In other 
words, the complex nature of wireless interference is not 
adequately captured by a conflict graph. On the other 
hand, the SIR-based interference model can be used to 
overcome this limitation. 

b. Channel model: It is assumed that the wireless channel 
is either to be static {i.e., not time-varying), or that 
the instantaneous CSI (Channel state information) at 
each time slot is available for scheduling collision free 
transmissions. However, wireless channels are seldom 
static due to fading, and the availability of CSI at each 
transmitter is not necessarily realistic in an adhoc setting. 

We present a brief summary of the assumptions made in the 
existing literature in the following table: 
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This work 

SIR 
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Stat. 
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« Inst. - Instantaneous channel gains are assumed to be 
known at each time slot. 

« Stat - Channel statistics (such as average channel gains 
or distribution) are assumed to be known. 

A time-varying channel is considered between the transmit¬ 
ter of a link and its corresponding receiver in fTl. However, 
the channel gains between the interfering links are assumed to 
be static. In 0, time varying channels are considered among 
all the links, and the interference is modelled by a conflict 
graph. However, the algorithm m can support only a fraction 
of the achievable rate region. A SIR model is considered in 
0 to a propose conservative algorithm that is suboptimal. An 
adaptive Aloha based algorithm is proposed in uni under time- 








varying channels. However the algorithm can only maximize Here \hii\'^R “ is the received power at the receiver in link i 
some utility functions and is not throughput optimal. from its intended transmitter and 


B. Our Contributions 

In this work, we consider a single-hop wireless network and 
propose a distributed scheduling algorithm. 

• We consider time-varying channels among all the links in 
the network. Further, the interference is modelled using 
the SIR model which is more realistic. 

• A key contribution of this paper is in the design of a 
Gibbs sampler Q based throughput optimal scheduling 
algorithm {Algorithm 1). In the algorithm we propose, 
each link only requires the average channel gains from 
its neighbouring links (defined later). In particular, in¬ 
stantaneous channel gains are not required, which makes 
our algorithm practical in a fast fading scenario, where 
the channel gains vary rapidly within a data slot. 

• We augment Algorithm 7, which allows only single 
link updates, and propose Algorithm 2, which performs 
parallel link updates and converges faster. 

The remainder of the paper is organised as follows. In 
Section m the network model is described. In Section Hm the 
spatial CSMA algorithm is presented and its throughput op¬ 
timality is proved. Numerical results are presented in Section 
EYl and we conclude in Section V. 

II. Network Model 

We consider a single-hop ad-hoc wireless network, and 
model the links using a bipole model introduced in ifTOl . In a 
bipole model, each transmitter is associated with a receiver that 
is at a distance R in some arbitrary direction. A transmitter and 
its corresponding receiver is referred to as a link. We assume 
that there are N links in the network. We use the set TV to 
denote all the links in the network. We assume that the time 
is slotted. 

We assume that the link distance R is much smaller than 
the distances of the transmitter and receiver to the other links. 
With this assumption, we can think of links as points in the 
Euclidean space (The results in this paper are not limited by 
this assumption. The assumption is taken to keep the expres¬ 
sions simple). Let rji denote the distance between the links 
i,j. We consider a standard path-loss model ||a;||““,a > 2. 
Channel model: The small-scale fading (power) between any 
pair of nodes is modeled by a unit power Rayleigh distribution 
and is assumed to be i.i.d across time and space. The channel 
gain between the transmitter of a link i and the receiver in link 
j is denoted by hji. Since hij is Rayleigh distributed, 
is exponentially distributed with unit mean. 

Interference model: A receiver successfully receives the packet 
of the corresponding transmitter if the received SIR is above a 
pre-determined threshold T. We consider interference limited 
networks, where the impact of thermal noise is negligible as 
compared to interference. Suppose M C J\f, be the set of 
links that are transmitting in the current slot. The SIR of a 
link i G A4 denoted by is given by. 


iGAt\{z} 

is the interference power from other concurrent transmissions. 
Queuing Dynamics: Each link has a separate arrival process 
and maintains its own buffer. denote the arrival rates 

of the links, [qi{t)]fLi denote the queue lengths of the links 
in time slot t. 

Assumptions on channel state information: We assume that 
each link knows the distances to its neighbours (defined later), 
the path loss exponent and the SIR threshold T. 

We now compute the probability that a transmission is 
successful in the presence of interference. 


A. Probability of successful link 

The probability of success for a link i G Ai denoted by 
Hi{M) is 


PiiM) = P(7i,Ai > T), 

= F {\hu\^ >R‘^TI{M\{i})) 
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where (a) and (c) follow from the exponential distribution 
of p and (b) follows from the independence of 

the fading variables. Let f(rij) := — , , c . Then the 

probability of success can be written as. 


n (1) 


Eor convenience, the probability of success is set to zero 
for the links that are not in the currently active set M., 
i.e., Pi[M) =0, Vz ^ M. Note that ([T]l is calculated, 
assuming all the active links in the network can contribute to 
the interference of a receiver. However, from the studies on 
statistical distribution of co-channel interference, the aggregate 
interference from the links beyond a certain distance can 
be safely neglected ED, ED- The radius beyond which the 
interference can be neglected is referred to as close-in radius, 
and is denoted by 77/. Hence for a link i, the interference 
from the active links outside a ball of radius 77/ around i can 
be neglected. Let A/i denote the set of links that are potential 
interferers of link i, i.e., the set of links within the ball of 
radius 77/ around link i. The links in Mi are referred to as the 
neighbours of link i. Thus, from E) the probability of success 
is. 


_ \hu\'^R /(Aj), Vie Ad, 

7(7W\{z}) j&Mi 
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where A4i := Mi H A4 is set of active links that are within 
the close-in radius of link i. Also note that, if none of the 
potential interferers of a link i are active, then it succeeds 
with probability one, i.e., 

= 1, if MnMi = 0. 

From (O, we can observe that the probability of success of 
a link depends only on the distances from its active neighbours 
Mli- This allows for computation of iJ,i{A4) by a simple 
neighbour discovery algorithm ini. 

We now characterize the capacity region in terms of the link 
success probabilities fj,i{Ml). 

B. Capacity Region 

Every subset of the links Ml C M, is associated with a TV¬ 
dimensional vector ti{Mi) = whose T-th element 

correspond to the probability of success of the link i (when 
Mi is the set of links that are transmitting). fJ,{Ml) can also 
be interpreted as the long-term rates that can be supported 
when the subset Ml is active. We refer to these vectors as rate 
vectors. 

The capacity region of the network is the set of all the arrival 
rate vectors for which there exists a scheduling algorithm that 
can stabilize the queues. It is known that the capacity region 
is given by 

A = {a e I 3e > 0, a(l + e) € Co (/r)}, 


• Neighbour discovery- Each link j G {z} U A/i executes a 
neighbour discovery 03 algorithm to compute the set of 
its active interferes in the previous slot, i.e., Mlj{t — 1). 

• Inactive weights- Each link j G Mi computes pj {Ml {t — 
1) \ {z}) from (|2l and subsequently computes the inactive 
weight 

:= g{qj{t)) fj.j{M{t-l)\{i}). (3) 

• Active weights- Link i obtains the inactive weights from 
its neighbours and computes the active weights as defined 
below. 


w] ■■= g {qi{t)) gt{Ml{t - 1) U {z}). 


• Update Probability- Link z computes its update probabil¬ 
ity p{t) as. 


p{t) = 


exp [wj) 


(4) 
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Link z chooses to transmit with probability p{t) and 
chooses not to transmit with probability 1 — p{t), i.e., 


M{t) = 


Mi{t — l)U{i} w.p. p{t), 
A4(f —l)\{z} w.p. l—p{t). 


Data slot: In the data slot, all the links j & Ml {t) will transmit. 


where, Co{p) is the convex hull of {p.{Mi)}MaN■ 

An arrival rate vector y G M" is said to be feasible if z/ G A. 
A scheduling algorithm is said to be throughput optimal, if 
the algorithm can stabilize the network for any feasible arrival 
rate. A maximum weight scheduling algorithm is known to be 
throughput optimal. In each time slot, the algorithm picks the 
schedule, Mi{t) = argmax P-j [Ml) qj{t). 

McN jeM 

Some of the notations used so far, are summarized below. 



Set of all the links in the network 

Ml{t)- 

Set of links that are active in slot t 

Mi- 

Set of potential interferers of link i. 

Ml,{t)- 

Set of active interferers of link i in slot t. 

Pi{Ml)- 

Rate of link i when the set of active 


links is Ml. 


III. Spatial CSMA 


In this Section, our distributed algorithm. Spatial CSMA is 
presented and its throughput optimality is proved. The key idea 
is to sample subsets (of links) so that sampled subsets provide 
a good approximation to the Maximum weight algorithm Q, 
©. Let g{x) be a real valued function of queue length. 
The details of the function g{x) are discussed subsequently. 


Algorithm!-. Spatial CSMA 


Intialization: Each link i G M pre-computes := f{rij) 
for all its neighbours j G Mi . 

Control slot: 

• Decision schedule- A link z G M, is picked uniformly at 
random. 


In Algorithm!, each time slot is divided into a control slot 
and a data slot. In the control slot, a link z is chosen at random 
(the implementation of this step is discussed later), and only 
this link is allowed to change its status (on/off) in this time 
slot. All other links will retain their status of the previous time 
slot. Link z and its neighbours execute a neighbour discovery 
algorithm to identify all their active neighbours. Lor example, 
the compressed neighbour discovery scheme 03 is a fast and 
efficient neighbour discovery algorithm which jointly detects 
all the active neighbours by allowing them to simultaneously 
report their identity. 

All the neighbours of the link z, use their neighbourhood 
information Mij{t — 1) computed in the previous step to 
calculate their rate vectors. Note that, all the links in Mlj{t—1) 
retain their status except for the possible change of the status 
for link z. Hence, to account for this possible change of the 
status of link z, we define two sets of weights namely inactive 
weights and active weights. The contribution of interference 
from link z is excluded for computing the inactive weights but 
included for computing the active weights. Link z uses these 
weights to compute its update probability p{t), and updates its 
status accordingly. In the data slot, all the active links transmit. 

A. Throughput Optimality 

Lemma 1. If the queue lengths are fixed at q = \qi]f^i, then 
Algorithm 1 corresponds to a Glauber dynamics Markov chain 
on the subsets Ml G M with a stationary distribution given 

by, 

n(A4) = ^ exp j ^ p.j (Mi) g(qj) j , \fMl cM, (5) 

XjeM J 







where Z is the normalizing constant. 

Proof: Proof can be found in Appendix- IVI-AI ■ 

If the queue lengths were indeed hxed (say at q) as required 
in Lemma [T] Algorithm! provides a good approximation for 
the maximum weight scheduler 0. This can be observed 
from as the stationary distribution 11 on the set of 
subsets, places the largest mass on the set Ai that maxi¬ 
mizes Pj which is precisely the max-weight 

j&M 

scheduler except for qj being replaced by g{qj). However, this 
replacement can be justihed if an appropriate function g{x) is 
chosen M- 

Lemma [T] assumes that the queue lengths are fixed. How¬ 
ever, the queue lengths are time-varying. Moreover, the time 
required for the Glauber dynamics to reach steady-state can be 
very long in general to assume that the queue lengths do not 
change. However, if appropriate slowly varying functions like 
log(0.1a::), loglog(a::-|-e) are used as g{x), it can be shown 0, 
m that Algorithm! does approximate the maximum weight 
scheduler in each time slot with a high probability and is hence 
throughput optimal. 

Lemma 2. If g{x) = log(0.1a::) or loglog(a:-|-e), the proposed 
spatial CSMA algorithm is throughput optimal. 

Proof: Follows from our Lemma [T] and Theorem 1, 
Proposition 2 in 0. ■ 

Remarks: While the techniques used are standard, the key 
contribution of this paper is the application of these techniques 
to design a throughput optimal scheduling algorithm for the 
SIR model with time-varying channels. 

Although Algorithm! is proposed for a Rayleigh fading 
model, it can be easily extended to other fading models without 
any additional effort. 

In 0, CSMA algorithm is considered on a conflict graph 
with time-varying link capacities. The authors of 0 show that 
the back-off parameter should have a exponential form of the 
channel gain. They obtain this by solving a maximum entropy 
problem. However, as we see from Lemma [T] the exponential 
form follows naturally from the max-weight formulation. 

B. Parallel Updates 

In the control slot of Algorithm!, it is assumed that a link 
(decision schedule) can be randomly picked at each time 
slot to update its status. However, the algorithm does not 
explicitly describe how to implement that step. Moreover, only 
one link is allowed to update its status in a given time slot. 
In Algorithm2, we relax the limitation of single-update, and 
provide a distributed algorithm to pick a decision schedule. 
The limitation of single-update can be relaxed by considering a 
block (parallel update) Gibbs sampler based algorithm, which 
allows for parallel updates and also converges faster. However, 
to ensure a distributed implementation of the (parallel update) 
Gibbs sampler, the set of links that can do parallel updates 
has to satisfy the following constraint. (See Lemma 0 for a 
formal proof.) 

If link i updates its status in a given slot, all the links whose 
current status information is being used in the computation 
of the update probability p(f) of link i, cannot update in the 
same slot. 


A set of links, which satisfy the above constraint is referred 
to as a decision schedule. The formal definition is as follows; 


Definition 1. Decision Schedule 

A set of links T> C M, is said to be a decision schedule if 


Mu 


U 


yjeN'i 


CM\V, 


Vi e V. (6) 


The intuition for this definition is as follows. From (|4|i, one 
can observe that the update probability p{t) of a link i, depends 
only on the weights of active links in M. Further, the weight 
of each link j € Afi, depends on the the status of its neighbours 
Mj. (See the computation of incative weights in Algorithml.) 
Hence, the required constraint translates to 0. 

Remark: This constraint is only on the set of links that can 
do parallel updates in a given slot. However, there are no hard 
constraints on the set of links than can transmit in a given slot. 
This is a key difference of this model compared to conflict 
graph based model. 

Generating a decision schedule V can be done in two steps. 
Step !: Generate a subset of links S, such that no two links 
in S are within the close-in radius of each other. 

Step 2: Initialize V to S, and update V by removing some 
links from V as follows. Each link k ^ S checks if any of its 
neighbours are in S. If more than one of its neighbours are 
present in S, then the neighbours are removed from V. 

The first step ensures M M\T>, while the second step 

ensures ( IJ Afj \ {* *} ) C M \ 22 so that 0 is satisfied. 

/ 

In 0 , a distributed algorithm is suggested for generating 
the subset S {step !). We extend that algorithm to generate 
the subset V. For the sake of completeness, we present the 
step ! from 0 . In Algorithm2, the control slot is divided into 
(W+2) control mini-slots for some W > 2. (This bound on W 
is to ensure that each link has a non zero probability of being 
selected in the decision schedule.) In each time slot, all the 
links in the network will execute Algorithm2 independently. 


Algorithm2: Decision Schedule Algorithm (at link i) 


Step!: Generating S 

• Link i selects a random (integer) backoff time Ti uni¬ 
formly in [0, W — 1] and waits for Ti control mini-slots. 

« If link i hears an INTENT message from a link in M 
before the {Ti + l)th control mini-slot, i will not be 
included in S and will not transmit an INTENT message. 

« If link i does not hear an INTENT message from any 
link in Mi, before the {Ti + l)th control mini-slot, it will 
send (broadcast) an INTENT message to all links in M 
at the beginning of the {Ti + l)th control mini-slot. 

- If there is a collision {i.e., if there is another link 
in M transmitting an INTENT message in the same 
mini-slot), link i will not be included in S. 

- If there is no collision, link i will be included in S. 

Step2: Generating V from S 

« If link i € S, it sends (broadcasts) an INTENT message 
to all links in M at the beginning of {W + l)th control 
mini-slot. 
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Fig. 1: Numerical Results 


• If link j ^ 5, it senses the channel for a possible collision 
{i.e., if more than one of its neighbours are in S, then 
all of them send INTENT messages which result in a 
collision) in {W + l)th mini-slot. 

- If there is a collision, link i will broadcast a DETECT 
message to all its neighbours in (W + 2)th control 
mini-slot. 

« If link i G 5 and it doesn’t hear a DETECT message in 
{W H- 2)th control mini-slot, it will be included in V. 


Lemma 3. If all the links in a decision schedule T) selected 
from Algorithm2, simultaneously update their schedules using 
Algorithml, then the stationary distribution of the resulting 
(parallel update) Glauber dynamics is given by (El). 

Proof: Proof can be found in Appendix IVI-BI ■ 

IV. Numerical Results 

In this Section, we evaluate the performance of the spatial 
CSMA algorithm. The results are compared to conflict graph 
based CSMA. This comparison requires the generation of an 
equivalent conflict graph for a given set of locations (of the 
links) as described below. 

Construction of conflict graph: In a conflict graph based 
interference model, each link in the network is represented 
by a vertex in a graph. Two vertices are connected by an 
edge, if their concurrent transmissions can possibly end up 
in a collision. Concurrent transmissions from any two links 
that are with in the close-in radius of each other can be 
unsuccessful (depending on the channel conditions). Hence, 
a pair of vertices are connected by an edge if they are in the 
close-in distance of each other. 

A. Simulation settings 

We consider a two dimensional square plane with side 
length 13. A homogeneous Poisson point process of density 
0.1 is generated. The generated points correspond to the 
locations of the transmitters. Each transmitter has its receiver 
at a distance of 0.25 in a random direction. The path loss 
exponent a is set to 2.5, the close-in radius Ri is set to 4, and 


the threshold SIR is set to 17 dB. The function g{x) (used in 
Algorithml) is set to log(0.1a;). In each time slot, the channel 
gains corresponding to unit power Rayleigh distribution are 
generated. 

B. Throughput performance 

In Eigure [T^ we illustrate the throughput performance of 
Algorithml, by plotting the average queue lengths to see 
for which arrival rates the system is stable. If the algorithm 
cannot stabilize the network for a given arrival rate, the queue 
length blows up. We consider homogeneous arrival rates for 
all the links. It can be observed that the SIR based algorithm 
supports a larger set of arrival rates compared to the graph 
based algorithm. This is because, in a conflict graph model, 
concurrent transmissions from two neighbouring links are 
strictly prohibited irrespective of the exact distance between 
them. However, in SIR model, the links make a better choice 
by considering the exact distances from its neighbouring links 
(thereby taking into account the severity of interference) while 
computing the update probabilities. 

C. Convergence rate 

In Eigure[T3 we compare the convergence rate of the spatial 
CSMA with the graph model by plotting the total queue 
evolution as a function of time. We consider a homogeneous 
arrival rate of 0.2 which is in the stable region of both these 
models. It can be observed that the queue reaches steady state 
much faster in the SIR model as compared to the conflict graph 
model. 

V. Concluding Remarks 

In this paper, we considered the SIR model with time- 
varying channels, and proposed a distributed CSMA algorithm. 
We further proved that the proposed algorithm is throughput 
optimal. We also proposed a parallel update algorithm with 
a better convergence rate. Using simulations, we observed 
that the SIR model supports a larger set of arrival rates, and 
converges much faster than the conflict graph based model. 


































VI. Appendix 

A. Proof of Lemma - 17 ] 

The Glauber dynamics (section 3.3.2 of ||4l) corresponding 
to the distribution 11, is a reversible Markov chain with state 
space {M \ Ai C Af}, and stationary distribution 11. The 
transition probabilities of that Markov chain are described 
here. From a given state Ai{t — 1), the chain moves to a 
new state as follows. A link i is chosen uniformly at random 
from JV and a new state is chosen according to the measure 11 
conditioned on the set of states in which status (on/off state) of 
all the links other than link i remain the same as in Ad (f— 1). 
In other words, the chain can only move to Ad(f — 1) U {*} 
or Ad (f — 1) \ {i} according to the conditional distribution 

n(Ad(f) I Ad(f) e {M{t - 1) U - 1) \ {*}}). (7) 

From (l5]l, it can be easily verified that the update probability 
p{t) in Algorithml correponds to the conditional distribution 
in Q. Hence, Algorithml corresponds to a Glauber dynamics 
Markov chain with stationary distribution given by 0. 

B. Proof of Lemma - 15] 

Let us represent a state Ad of the Glauber dynamics Markov 
chain with a N— dimensional binary vector a whose elements 
are given by ai = \{i G Ad}, i G M. Thus, the equivalent 
state space of the Markov chain is {u | u € {0, 1}^}. a{t) 
can be shown as a Random field (chapter 7 in 0) on M 
with {ai{t)}i^jif being the underlying random variables. We 
construct an undirected graph G(y, E) with V = N and edges 
given by the following definition of the neighbourhood on G. 



where, Na{i) is the set of neighbouring vertices of vertex 
i in the graph G. Let us denote, a{S) = [iTiJigs for 
S C A/". Observe that the conditional distribution, ai (t) given 
a{JV \ {i})(f) corresponds to the distribution in 0 (which 
corresponds to the update probability p{t) as discussed in the 
proof of lemma 0. Hence, 

V{a^{t) = l\a{^^\{i}){t))=p{t). ( 8 ) 

From one can observe that the update probability p{t), 
depends only on the weights of active links in A/}. Further, 
the weight of each link k G Ni given by 0, depends on the 
the status of its neighbours A/}- In other words, the update 
probability pit) of link i, depends only on the status of the 
links in No{i)- Using this observation in ® gives, 

P(cri(f) I o-(A/'\ {z})(f)) =¥{aiit) I cr(AG(*))(f)) ■ (9) 

In other words, the random variable aiit) is independent 
of all other random variables given the random variables 
corresponding to its neighbours in G, i.e., ait) is a Markov 
random field 0 with respect to G. 

In a parallel update Glauber dynamics corresponding to H, 
if T> is the set of links that are selected for update, then the 
update rule is specified by the following joint law 0, 

FiaiV)it) I aiM \ V)it)) = n(a(f) | a(Af \ 77)(f)). (10) 


From 0, i gV ^ Ncii) C A/"\X>. Using this property along 
with 0, we can write (fTOl i as a product of single-site update 
probabilities as follows. 

n^iV)it)\aiM\V)it))= n P(a.(f)|a(ArG(*))(f)). 

iGV(t) 

In other words, we can use the same update rule as in 
Algorithml for all the links in X>(f) simultaneously and still 
converge to the required stationary distribution given by 0. 
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